THE INTEGRALS IN GRADSHTEYN AND RYZHIK. 
PART 8: COMBINATIONS OF POWERS, EXPONENTIALS AND 

LOGARITHMS. 
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Abstract. We describe some examples of integrals from the table of Grad- 
shteyn and Ryzhik where the integrand is a combination of powers, exponen- 
tials and logarithms. The expressions for some of these integrals involve the 
Stirling numbers of the first kind. 



1. Introduction 

The uninitiated reader of the table of integrals by I. S. Gradshteyn and I. M. 
Ryzhik |4J will surely be puzzled by choice of integrands. In this note we provide 
an elementary proof of the evaluation 4.353.3 



(I) [\ax + n + l)x n e ax \nxdx = e a Y j {-l) k - 1 - + " ! 

Jo tZn {n-k)la k+1 



n+l ' 

k=0 



We also consider the integrals 

(2) q n := / x n e 

Jo 

and the companion family 

(3) Pn := f x n e 

Jo 



The integral q n corresponds to the case a = —1 in (fT]). Section [3] provides closed- 
form expressions for p n and q„ . Section [4] considers the generalization 



(4) Pn{a) 



f x n e- ax dx and Q n {a) = f 
Jo Jo 

The main result of this section is the closed-form expressions 
(5) P n (a) := / x n e— dx = — T 1 - e~ a £ 



o 



k=0 



and 

fc-l 

/ ■ n< , — . i i 

Qn(a) ■= 



r 1 i n i / fc-i j \ 

^0 j — i \ „-_n J' I 



fe=l \ j=0 
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where 

r 1 i 

(6) Q (a)=/ e- ax lnxdx = — (7 + In + r(0, a)) , 

Jo a 

and r(0, a) is the incomplete gamma function defined by 

/>oo 

(7) Y{a,x) := / t^e'* dt. 

2. The evaluation of 4.353.3 

The identity 

(1) 4- (x n+1 e ax ) = (ax + n + l)x n e ax 
ax 

and integration by parts yield 

(2) f {ax + n + l)x n e ax \nxdx = - [ x n e ax dx. 
Jo Jo 

This last integral appears as 3.351.1 in 0]. We have obtained a closed-form ex- 
pression for it in [5] • A new proof is presented in Section 3J 

A closed form expression for the right hand side of ([2]) is obtained from 

1 < d 



(3) / x n e ax dx , , 

ota 

The symbolic evaluation of © for small values of n E N suggests the existence of 
a polynomial P n (a) such that 

The next lemma confirms the existence of this polynomial. 
Lemma 2.1. The function P n (a) defined by 

is a polynomial of degree n. 

Proof. Let D = 4-. Then D n+1 = D(D n ) produces the recurrence 

(6) P„+i(a) = aP;» + (a - n - l)P„(a). 

The initial condition Pa(a) = 1 and © show that P n is a polynomial of degree 
n. □ 

Theorem 2.2. The polynomial 

(7) Q n (a) := (-l) n P n (-o) 
has positive integer coefficients, written as 

n 

(8) Q n (a) -^&„, fc a fc . 

k=0 
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These coefficients satisfy 

(9) K+i,o = {n + 1)671,0 

6 n +i,fe = (n + l-fe)6„,fc + 6 nife _i, l<fc<n 

Moreover, the polynomial Q n (a) is given by 

n 

(10) 

Proof. The recurrence (jSJ) yields 

(11) Q„+i(a) = -aQ' n {a) + (a + n+ l)Q n (a). 

The recursion for the coefficients b n: k follows directly from here. Morover, it is clear 
that b ntn = 1 and b n> o — n\. A little experimentation suggets that b n _k = n!/fc!, 
and this can be established from ([9]). □ 

This proposition amounts to the evaluation of 3.351.1 in 0]: 

(_-\\n+l „| p au ™ I 



(12) x n e ™ dx = ^l_Ifi + « V^(-l) 

/„ a n+l a n+l fcl ; 

fc=0 



u a 







The reader will find a proof of this formula in [2]. 

3. A NEW FAMILY OF INTEGRALS 

In this section we consider the family of integrals 
(1) q n := / x n e~ x \nxdx, 



and its companion 

(2) 







Lemma 3.1. The integrals p n , q n satisfy the recursion 

(3) p n+1 = {n+l)p n -e~ x 

(4) q n +\ = (n+l)q n +p n 

Proof. Integrate by parts. □ 
The initial conditions are 

(5) po — 1 — e _1 and q Q = / e~ x In xdx = 7 — Ei(— 1). 

Jo 

Here 7 is Euler's constant defined by 

™ 1 

(6) 7 := lim > - — lnt 
with integral representation 

(7) 1=1 er x \nxdx 



n — >oo L — ' \z 

k=l 



n 
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given as 4.331.1. The reader will find in [3] a proof of this identity. The second 
term in ([5]) is converted into 



e 



(8) / e x In x dx = / dx 

Ji Ji x 

and this last form is identified as Ei(— 1), where Ei is the exponential integral defined 

by 



(9) Ei(z) = -^ 



dx. 



In the current context, the value of Ei(— 1) will be simply one of the terms in the 
initial condition q . 

We determine first an explicit expression for p n . The recursion §3§ shows the 
existence of integers a n , b n such that 

(10) Pn=a n + b n e~ 1 , 
with ao = 1, bo = — 1. From ((3|) we obtain 

(11) a n+1 + b^e- 1 = (n + lR + [(n + 1)6„ - 1] e -1 . 
The irrationality of e produce the system 

(12) a n +i = (n+ l)a n , with a = 1, 

(13) &n+i = (w + l)&n — Ij with bo = — 1. 

The expression a„ = n! follows directly from (fT2|) . To solve (fT3|) . define i?„ := 6„/n! 
and observe that 

1 



(14) B n+1 =B n - 

[n + 1) ! 

that telescopes to 

n 1 

(15) 6 « = -" ! Efc|- 

fe=0 

We have shown: 

Proposition 3.2. The integral p n in (0) is given by 

(16) ^f^^T^i)' 

We now determine a similar closed-form for g„. The recursion (j4|) shows the 
existence of integers c„, d n , fn such that 

(17) q» = c„ + d„e _1 + / n g - 

In order to produce a system similar to (|12I13[) we will assume that the constants 
1, e" 1 and qo = —(7 + Ei(— 1)) are linearly independent over Q. Under this as- 
sumption (j4|) produces 

(18) Cn+l = (n+l)c n + n!, 

(19) d„+i = (n + l)c ra - niy^ — , 

fc=0 

(20) / ti+1 = (n +!)/„, 
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with the initial conditions Co — 0, do = and /o = 1. 

The expression f n = n>\ follows directly from (f20l) . To solve (fl8|) and (fT9|) we 
employ the following result established in [T] . 

Lemma 3.3. Let a n , b n and r n be sequences with a n , b n ^ 0. Assume that z n 
satisfies 

(21) a n z n = b n z n ^i + r n , n > 1 

with initial condition Zq. Then 



(22) z n 

We conclude that 
(23) 

and 
(24) 



bib 2 ■ ■ - b n 
a\a 2 ■■■an 



E 

k=l 



a\02 ■ ■ ■ ak-i 
bib 2 



-ru, 



n 1 

'V- 



fe=l 



n A; — 1 



fc=l j=0 



The expression for c n shows that they coincide with the Stirling numbers of the 
first kind: c n = \s(n + 1, 2)|. 

We have established 
Proposition 3.4. The integral q n in fip is given by 



(25) 



x n e x \nxdx — n\ 



1 ™ 1 

-T- 

e ^ k 

fc=i 



k-l ^ 

J=0 



<lo 



Example 3.5. The expressions for p n and q n provide the evaluation of 4.351.1 in 

(26) / (1 ~ x)e- x lnxdx 



1-e 

e 

by identifying the integral as go — 1\ ■ The recurrence (HJ) shows that 
(27) q - qi = -p = e^ 1 - 1, 

as claimed. 

Example 3.6. The evaluation of 4.362.1 in 0] 



(28) 



xe x - x) dx = / (1 - i)e 1-i lntdt 



is achieved by observing that this integral is e(qo — q±) = 1 — e. 
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4. A PARAMETRIC FAMILY 

In this section we consider the evaluation of 

(1) Pn{a) := 

(2) Q n (a) := 

The integrals q n considered in Section[3]corresponds to the special case: q n — Q n (l). 
We now establish a recursion for Q n by differentiating ((2]). 

Lemma 4.1. The integral Q n (a) satisfies the relation 

(3) Q n+ x(a) = -4-Qn(a). 

da 

To obtain a closed- form expression for Q„(a) we need to determine the initial 
condition 

(4) Q (a) = / e- ax \nxdx. 



This is expressed in terms of the incomplete gamma function defined in 8.350.1 by 

/>oo 

(5) T(a,x) := / i^Vdt. 



Observe that r(a,0) = T(a), the usual gamma function. 
Lemma 4.2. The initial condition Qo(a) is given by 

f 1 1 

(6) Qo(a)= e- ax lnxdx= — (7 + lno + T(0,a)) . 

Jo a 
Proof. The change of variables t — ax yields 

(7) Qo(o) = - / e- t \ntdt- — (l-e- a ). 

a J a 

Then 



(8) J e^lntdt^ / e^lntdt- / e^hxtdt. 
The first integral is 

(9) / e^hxtdt 



that simply reflects the fact that 7 — — r'(l). Integrating by parts yields 



(10) 



/ e'Hntdt = e- a lna + T(0,a). 

J a 



The formula for Qo(a) is established. □ 

We now determine a closed- form expression for P n (a) and Q n (a) following the 
procedure employed in Section [3] 
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Lemma 4.3. The integrals P n and Q n {a) satisfy the recursion 



(11) 
(12) 



1 



Pn+iO) = - ((n+ l)P„(a) - e- a ) 
a 

n+iO) = -((« + l)Qn(a) + P„(a)) . 



( 7 + r(0,a) +lna). 



TTie initial conditions are given by 

(13) P (a) = -(1 - e~ a ), and Q (a) = — 

a a 

Proof. Integrate by parts. 

We conclude that we can write 

(14) P n (a) = A n (a) - B n (a)e- a 7 
and 

(15) Q n (o) = C„(a) - £> n (a)e- a - £„(a)( 7 + T(0, a) + In a). 
Lemma 4.4. TTie recursions and M2\) imply that 



a 



(16) A n+1 (a 
B n +i(a 

Cn+l( a 

D n+1 (a 

with initial conditions 

(17) A (a) = B {a 



(n + l)A n (a), 



= -[(n + l)B„(o) + l], 
a 

= - [(n + l)C„(a) + A„(a)], 
a 

= i[(n + l)D„(a)+fl B (o)], 
a 

= -(n+l)S n (a) 
a 



So (a) 



and Co(a) = -Do(a) = 0. 



These recursion can now be solved as in Section [3] to produce a closed-form 
expression for the integrals P n (a) and Q n (a). We employ the notation 



(18) 



n 



fe=i 



for the harmonic numbers and 

n 

(19) Ex Pn (x)=Y / 



k=Q 



X 

71 



for the partial sums of the exponential function. 
Theorem 4.5. Let a e R one? neN. TTien 



(20) 
and 



P n {a) := J x n e~ ax dx = [l - e- a Ex Pn (a)} , 



Qn(a) 



In a; dx 



,n+l 



fe=i 
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where G(a) = —aQo(a) — 7 + r(0, a) + In a. 

These expressions provide the evaluations of two integrals in [4]. 
Example 4.6. Formula 4.351.2 states that 

r 1 1 

(21) / e- ax (-ax 2 + 2x)\nxdx = [-1 + (1 + a)e" Q l . 
Jo a 

In order to verify this, observe that the stated integral is 

(22) -a [ x 2 e~ ax \nxdx + 2 [ xe~ ax \nx dx = -aQ 2 (a) + 2Q 1 (a). 

Jo Jo 

The expressions in Theorem 14.51 now complete the evaluation. 
Example 4.7. Formula 4.353.3 in [Jj gives the value of 

(23) I n (a):= f (-ax + n + l)x n e - ax In a: dx. 

Jo 

Observe that 

(24) I n (a) = -aQ n+l (a) + [n + l)Q n (a), 

and using the recursion (|12[) we conclude that I n (a) = —P n (a). The expression in 
Theorem 14.51 is precisely what appears in [1] . 

We conclude with the evaluation of a series shown to us by Tewodros Amdeber- 
han. Expand the exponential term in (|21|) and integrate term by term to obtain 
(25) 

x f \ k J / 7i 1 \ 

£ fc! (n+iW = H n + 1} + 111 a + m a) + era E ^ ExPfc - i(a) • 

fc=0 \ ' ' I \ fc =0 / 

Here 

(26) 1>(x) = £M 

T(x) 

is the digamma function defined in 8.360.1 of 4J. the identity 

(27) ^(n+l) = fln- 7) 

that is a direct consequence of the functional equation T(x + 1) = xT(x) and 
r'(l) = —7, was used to transform (|25|) . 



The identity (|25| can be used to provide multiple expressions for the incomplete 
gamma function, such as 

(28) r ,t ; )>n T" + ^(n + 1) - m. - e- 1 a y > , 

y ' J a x ^ nlkl(n + l + k) 2 y ' ^ k 

k—0 k=l 

and the special case for n = 0: 

(_l)fc fc+i 



(29) / dx = -7-lna + ^ 



x ^ (fc + l)!(fc + l) 



These issues will be explored in a future publication. 
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